Concentration Dependence of the Effective Mass of 3 He Atoms in 3 He- 4 He Mixtures 
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Recent measurements by Yorozu et al. (S. Yorozu, H. Fukuyama, and H. Ishimoto, Phys. Rev. 
B 48, 9660 (1993)) as well as by Simons and Mueller (R. Simons and R. M. Mueller, Czhechoslowak 
Journal of Physics Suppl. 46, 201 (1976)) have determined the effective mass of 3 He atoms in a 
3 He- 4 He-mixture with great accuracy. We here report theoretical calculations for the dependence 
of that effective mass on the 3 He concentration. Using correlated basis functions perturbation 
theory to infinite order to compute effective interactions in the appropriate channels, we obtain 
good agreement between theory and experiment. 



The effective mass of a 3 He atom in liquid 4 He is due to 
two effects. The dominant one is hydrodynamic backflow 
of the 4 He liquid around the impurity [0J|]. This effect 
is largely independent on the 3 He concentration. The 
second effect is the dynamics imposed on the 3 He com- 
ponent by the Pauli principle ||. This effect causes a 
noticeable concentration dependence which has recently 
been measured with high accuracy by Yorozu et al. Q 
and by Simons et al. it is the subject of this paper. 

Ground state properties of 3 He- 4 He mixtures like the 
energetics of the system and its local structure are today 
quite well understood from a microscopic point of view 
H . For a microscopic understanding of the hydrodynamic 
effective mass, rather advanced methods are needed Q] 
due to the high density of the 4 He background. On the 
other hand, the 3 He component in the mixture is dilute, 
and the interaction between individual 3 He atoms is dom- 
inated by phonon exchange. Therefore, much simpler 
methods than those necessary || for obtaining reliable 
results in pure 3 He are adequate for calculating ground 
state properties of dilute mixtures || . 

Microscopic many-body theory postulates an empirical 
Hamiltonian 



H 
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a/3 i,j 



(1) 

that contains only a two-body interaction || and the 
masses of the two types of particles. One then uses the 
Feenberg form |l(],[ll| for the ground state wave function, 

* ({r^})=e^^ Q, »<i>o({rf)}), 

N n ,N„ 



a/3 i,j 



(2) 



afil i,j,k 



where $o({ r i }) is a Slater determinant of plane waves 
ensuring the antisymmetry of the Fermion component 



of the mixture. The superscripts a, /?, ... refer to the 
type of correlated particles; the prime on the summa- 
tion symbol in Eqs. (|l|), (||) indicates that no two pairs 
(i,ct), (J,/3) can be the same. The correlation func- 
tions u a,3 {Ti,Tj) and u Q/37 (rj, rj, rfe) are determined by 
the functional minimization of the ground state energy 



E a = (* |ff|*o>/(*o|*c 



(3) 



with respect to the correlation functions |13-H . Details 
of the procedure have been described in Ref . [p| . 

Within the above variational theory, the single par- 
ticle excitation spectrum is calculated by allowing for 
an occupation of single particle orbitals in the Slater 
function that is different from the Fermion ground state 
nl = 9(kp — k), where kp is the Fermi momentum of 
the 3 He component. The single particle spectrum can 
then be calculated by the variation 



^{k)= 5 -^ 



t 3 (k)+u(k) + U 



(4) 



where we use t a (k) = h 2 k 2 /2m a for the kinetic energy 
of a free particle of species a. Uo is a constant related 
to the chemical potential, and u(k) a momentum depen- 
dent average field. The general graphical analysis of the 
variational single-particle field u(k) has been carried out 
in Ref. fL5|| ; if the 3 He component is dilute, the average 
field can be written in the form of a Hartree-Fock field 
in terms of a local, effective interaction W e s(q) 



u{k) 



d 3 q 



(2^) 3 p 3 q 



n<°> Weff(|k-q|) 



(5) 



It is known |q] that the nai've use of a variational wave 
function of the type (^) leads to an effective mass of pure 



3 He that is — in sharp contrast to experiments |L6 17 
— less than one p8fl. The cause of this deficiency is that 
the wave function (0)) describes the average correlations 
between particles, but it is not sensitive to the specifics 
of the correlations in the vicinity of the Fermi surface. 
The cure for the problem is Correlated-Basis Functions 
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(CBF) theory g(J to infinite order |9[. The theory can 
be mapped on a Green's function approach in terms of 
effective interactions p(| that are provided by the varia- 
tional theory. The single particle properties are described 
by a complex self-energy £(£;, u>); the single particle spec- 
trum e(k) is obtained from the solution of equation 



e(k) = t 3 (k) + 'E{k,s(k)) 



(6) 



If only one-phonon coupling processes are considered, 
the self-energy T,(k,E) is given by the so-called GOW- 
approximation |2lLp2| 



d 3 qd(huj 
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G (0) (fc,w) = 
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huj — t 3 (k) + in Tilj — t 3 {k) — in 
is the free single-particle Green's function and 

V eS (q,u) = Vp 3 -h(<z) + E %(?)Xa/j(9,«)%(?) ( 9 ) 

a/3 

is the effective, energy dependent 3 He- 3 He interaction. In 
Eq. (|9|), Vpf h (q) is the local, particle-hole irreducible in- 
teraction matrix ||, and X<*,8(?> k-O is the density-density 
response matrix. The particle-hole irreducible interac- 
tions, which are in conventional Green's functions theo- 
ries the most significant source of uncertainty, are pro- 
vided by the variational ground state theory. 

To separate the "hydrodynamic" and the "fermionic" 
component of the self-energy, we rewrite the single- 
particle Green's function as 



G (0) (fc,w) 



1 



hcj — t 3 (k) + in 
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ts{k) — fiw ~ if] — h*- + if) 

(10) 



= Gg»(* |W ) + G™(*,w) 



and, correspondingly, the self-energy as 
£(&, E) = Sh(A, E) + S F (fc, £■) 



(11) 



The "hydrodynamic" part E 1 ) of the self-energy is, 

in the limit p 3 — > 0, identical to the self-energy of a single 
3 He impurity in the 4 He host liquid (c/. Eq. (3.30) of 
Rcf. 0); since we have found in Ref. |(| that the con- 
centration dependence is generally weak it is appropriate 
to identify T,n(k,E) with the self-energy of a single 3 He 
atom also at finite concentrations; we will return to this 
quantity later. 

Let us first focus on the second term in Eq. (pi]). The 
energy integration yields the compact form 



J 0-n^V eS O,- q ,E-t 3 ( q )) 



(12) 



The relationship to Eq. (@) is apparent: If V c s{q,io) were 
energy independent, Yip (k, u>) would be just a Hartree- 
Fock average field of the form (||). Indeed, the expres- 
sion (JsJ) can be derived from Eq. ( [i"2"|) using the same 
"average-energy" procedure that has been employed to 
establish the connection between the parquet-diagram 
theory and the optimized HNC theory p3|, namely to 
identify 



W eS {q) = V eS {q,Q(q)) 



(13) 



where the average energy hQ(q) is chosen such that both 
the energy dependent and the energy independent inter- 
actions produce the same static structure function. Wc 
stress that this result is only an observation on how the 
static approximation and the G0W approximation are re- 
lated. It does not imply that this approximation is also 
adequate for the single-particle properties. 

With the relationship between theories established, we 
now turn to the numerical application. The ingredients 
of the theory — the effective interactions V™_ h (q) and 
W e s(q) as well as the Feynman spectrum e^ 4 '(k) have 
been obtained in Ref. The first important quantity 
is the hydrodynamic effective mass. Due to the high den- 
sity of the background more elaborate methods than the 
G0W approximation must be used for a quantitative pre- 
diction Q. To identify the Fermi-liquid effects we are 
interested in here, we have, however, not used these re- 
sults, but rather let the hydrodynamic mass be a free 
parameter. After the concentration dependence was cal- 
culated from the Fermi-liquid contributions we made a 
single parameter fit to the experiments of Refs. jl and || 
to optimize the overall agreement. That way, we arrived 
at the following interpolation formulas for the hydrody- 



namic mass: 



m ?> 



— = 2.171 + 2.334r + 4.547r^- 9.329r 3 (14) 



cxpt 



for the data of Ref. fl4fl and 



m 3 / cxpt 



2.124+ 2.630r + 2.000r 2 



(15) 



from those of Ref. ||. Here, r = pi/po — 1, pi is the 4 He 
density and /3o=0.02183A -3 is its value at the satura- 
tion vapor pressure. Typically, the discrepancy between 
the two different extrapolations is 0.03, these values are 
throughout the full density regime about 0.1 above our 
theoretical calculation of Ref. [J7J . 

A possible complication to be considered is the mo- 
mentum dependence of the hydrodynamic effective mass 
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since the particles at the Fermi surface have a finite mo- 
mentum. For that purpose, we have calculated the single- 
impurity spectrum, 



hu{k) = t z {k) + Y, H {k 1 uj{k)) 



(16) 



and have determined the momentum dependence of the 
hydrodynamic mass by writing, in a momentum regime 
< k < 0.4 A (chosen such that the upper bound- 
ary is comparable to the Fermi momentum at the highest 
pressure and concentration) the spectrum in the form 



ZrriH 



bk l 



(17) 



We have verified that this form is adequate at small 
momenta and has a weak density dependence in agree- 
ment with experiments. |^J,^5| We have used the value 
b = —0.074 at all densities. 

Three calculations have been carried out to determine 
the Fermi-liquid contributions to the effective mass of 
the 3 He component as a function of concentration and 
density. The transformation from the density and con- 
centration dependence to the pressure and concentration 
is done using experimental results given in Ref. p6| . 
The first is the simple approximation (^|) . To account for 
the hydrodynamic backflow, one must supplement the 
Fermion contribution (^) by the hydrodynamic contribu- 
tion; the spectrum has the form 

e {3) (k) = hcj H (k)+u(k) + U (18) 

where the Fermi correction u(k) is given in Eq. (|s|). The 
effective masses derived from this spectrum are signifi- 
cantly above the experimental one, cf. Fig. |l|. 

In the next step, we use the full self-energy in an "on- 
shell" approximation 



4to 



(19) 



This form of the self-energy relaxes the approximations 
made by the variational theory, we see in Fig. [j] that the 
agreement with the experiment is no better than in the 
FHNC approximation; the on-shell approximation, the 
effective mass now being significantly below the experi- 
mental value. 

Finally, we carry out a self-consistent calculation of the 
effective mass. Due to the low concentrations, we may 
assume a single-particle spectrum of the form ts(k) = 
fiujH(k) in the Green's function (j^) and, consequently, in 
Eq. (12); note that the hydrodynamic mass is included 
in the Green's function. This effective mass is then de- 
termined self-consistently by requiring that the spectrum 
e^(k) determined by 



,(3) 



(k) = hu>H(k) + Hp (k, Tvjju{k)) 



(20) 



can be fitted by the same effective mass that has 
been used in the self-energy. This theoretical calcula- 
tion produces, especially at lower densities, a slightly 



stronger concentration dependence than seen experimen- 
tally. Compared with both the FHNC and the "on-shell" 
approximation, the self-consistent result shown in Fig. |l| 
appear quite satisfactory. 

The essential difference between the on-shell and the 
self-consistent calculation is that information about the 
hydrodynamic backflow has gone into G { p\k,u;). Note 
that hydrodynamic backflow is — as a feature of exci- 
tations — generically not included in the approximation 
(||). Indeed, generalizing the variational theory to dy- 
namic correlations without introducing the coupling to 
the background would lead to an effective mass close to 
the on-shell mass, which is notably too low. 

To produce Fig. [l] we have used - as stated before 
- the hydrodynamic mass given in Eqs. (|l4|) or (jlH), 
respectively. Our calculations predict a curvature of the 
effective mass as a function of concentration, leading to a 
hydrodynamic effective mass slightly lower than the one 
predicted by linear extrapolation, cf. Table |[ Such a 
curvature is caused by the Fermi- functions, already the 
simple approximation (|l^) predicts a behavior 

m*(x) = ran + ax 2/3 + bx + cx 5/3 + dx 7/3 . (21) 

The precise values of these coefficients is a matter of a 
microscopic calculation. In Table | we list the values of 
a, 6, and c for different pressures as obtained from the 
least square fit to the self-consistent solution of Eq. fl2C|). 

We have demonstrated that the concentration depen- 
dence of the effective mass is essentially a Fermi-liquid 
effect, enhanced by the renormalization of the single par- 
ticle propagator through hydrodynamic backflow. The 
technical simplifications caused by the relatively low den- 
sity of the Fermion component has allowed to highlight 
the relevant physical mechanisms quite clearly. The data 
of Ref. H are consistently below those of Ref. ||; the 
difference of 0.03 might be attributable to the different 
pressure gauge. This difference might appear negligible, 
but we need to point out that it produces uncomfortably 
large uncertainties in predictions for the first antisym- 
metric Landau parameter Fq from magnetic susceptibil- 
ity measurements J27]] . 

On the other hand, we find that the concentration 
dependence of the effective mass — in other words the 
generic Fermi-liquid effect — is consistent between both 
sets of experiments and theory. Extrapolations to zero 
concentration also appear to be consistent. 

Besides providing an accurate microscopic calculation 
of the concentration dependence of the 3 He effective mass 
in He, we have analyzed various procedures for calcu- 
lating the effective mass. By comparing the static calcu- 
lation with the dynamic CBF calculation, we have dis- 
cussed an instructive example for the delicate interplay 
between single-particle and Fermi-liquid effects in 3 He- 
4 He mixtures. 
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FIG. 1. Theoretical and experimental effective mass ratio 
m* (P, x) /m as a function of pressure P and concentration 
x. The thick full curve is the fully self-consistent result with 
the hydrodynamics mass fitted to the experiments of Ref. Q] 
(circles with error bars). The thin, dashed-dotted curve is 
the on-shell approximation and the thin, short dashed curve 
the static approximation using the same hydrodynamic mass. 
The thick dashed curve is the fully self-consistent result using 
the hydrodynamic mass fitted to the results of Ref. D (boxes 
with error bars). 



TABLE I. Pressure dependence of the hydrodynamic effective mass and the coefficients of the expansion (Ell) as obtained 
from the present calculation and fitted to the data of Ref. Q. Also shown are the linear extrapolation of Ref. [^as well as our 
fit to these data. 



P (Atm) 


tuh (this work) 


a 


b 


c 


d 


run Ref. [ 




rriH fitted to Ref. 





2.18 


1.49 


1.39 


- 18.2 


36.7 


2.23T0.02 


2.15 
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2.31 


1.07 


3.00 


- 22.6 


40.2 
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2.44 


0.789 


4.48 


- 28.2 


50.4 


2.52T0.02 


2.39 


15 


2.54 


0.501 


6.17 


- 36.1 


66.8 








20 


2.64 


0.310 


7.41 


- 42.1 


80.1 


2.70T0.03 


2.62 
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